Astronomy & Astrophysics manuscript no. 021004-ac 


©ESO 2010 


April 8, 2010 





GRB 021004: Tomography of a gamma-ray burst progenitor and 

its liost galaxy * 

- - - A. J. Castro-Tirado\ P. M0ller^, G. Garcia-Segura^ , J. Gorosabel^ E. Perez\ A. de Ugarte Postigo^, E. Solano^, 
2 D. Barrado-Navascues^, S. Klose^, D. A. Kann^, J. M. Castro Ceron'', C. Kouveliotou*^ , J. P. U. Fynbo^, J. Hjorth^, 
^ ; H. Pedersen^ E. Pian^-"'-", E. Rol^^.n^ g paiazzi^'^, N. Masettii"*, N. R. Tanvir^^ P M. Vreeswijk^ M. I. 

. Andersen^, A. S. Fruchter^^, J. Greiner^^, R. A. M. J. Wijers^^, and E. P. J. van den Heuvel'^ 

U '. 

' Institute de Astrofisica de Andalucfa (lAA-CSIC), Glorieta de la Astronomia s/n, E- 18.008 Granada, Spain. 
- European Southern Observatory, Karl-Schwarzschild-Strae 2, 85748, Garching bei Munchen, Germany. 
Institute de Astronomia, Universidad Nacional Autonoma de Mexico, Apdo. Postal 877, Ensenada 22800, Baja California, Mexico. 
INAF, Osservatorio Astronomico di Brera, via E. Bianchi 46, 23807 Merate (LC), Italy. 
I I , ' Laboratorio de Astrofisica Estelar y Exoplanetas, Dpto. Astrofisica, Centro de Astrobiologia (CSIC/INTA), P.O. Box 78, 28691 Villanueva 

de la Caiiada (Madrid), Spain. 
^T^ ' Thiiringer Landesstemwarte Tautenburg, Sternwarte 5, D-07778 Tautenburg, Germany. 

European Space Agency (ESA), European Space Astronomy Centre (ESAC), P.O. Box - Apdo. de correos 78, 28691 Villanueva de la 
Canada, Madrid, Spain. 



1^ 



^ ** NASA Marshall Space Flight Center, NSSTC, 320 Sparkman Drive, Huntsville, Alabama 35805, USA. 

^ ' Dark Cosmology Centre, Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, DK-2100 K0benhavn 0, Denmark. 

^ INAF - Osservatorio Astronomico di Trieste, via Tiepolo, 11, 34131 Trieste, Italy. 

C/3 . " Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, 1-56126 Pisa, Italy. 

Department of Physics & Astronomy, University of Leicester, University Road, Leicester, LEI 7RH, UK. 

'^^ Astronomical Institute 'Anton Pannekoek", University of Amsterdam, PO number 94249, 1090 GE, Amsterdam, NL. 

^ INAF - lASF di Bologna, via Gobetti 101, 1-40129 Bologna, Italy. 

>- ' '5 Space Telescope Science Institute, 3700 San Martin Dr, Baltimore, MD 21218-2463, USA. 

T-H ' Max-Planck-Institut fiir extraterrestrische Physik, 85748 Garching, Germany. 

m ' 



^ ^ Received 24 December 2009; accepted 26 March 2010 

■ ABSTRACT 

o : 

, Aims. Analyse the distribution of matter around the progenitor star of gamma-ray burst GRB 021004 as well as the properties of its host galaxy 
with high-resolution echelle as well as near-infrared spectroscopy. 
^ Methods. Observations were taken by the 8.2m Very Large Telescope with the Ultraviolet and Visual Echelle spectrograph (UVES) and the 
• Infrared Spectrometer And Array Camera (ISAAC) between 10 and 14 hours after the onset of the event. 

Itn] Results. We report the first detection of emission lines from a GRB host galaxy in the near-infrared, detecting Ho- and the [O III] doublet. These 
allow an independent measurement of the systemic redshift (z = 2.3304 ± 0.0005) which is not contaminated by absorption as the Lyo- line 
is, and the deduction of properties of the host galaxy. From the visual echelle spectroscopy, we find several absorption line groups spanning a 
range of about 3,000 km s"' in velocity relative to the redshift of the host galaxy. The absorption profiles are very complex with both velocity- 
broadened components extending over several 100 km s"' and narrow lines with velocity widths of only ~ 20 km s"'. By analogy with QSO 
absorption line studies, the relative velocities, widths, and degrees of ionization of the lines ("line-locking", "ionization-velocity correlation") 
show that the progenitor had both an extremely strong radiation field and several distinct mass loss phases (winds). 

Conclusions. These results are consistent with GRB progenitors being massive stars, such as Luminous Blue Variables (LBVs) or Wolf-Rayet 
stars, providing a detailed picture of the spatial and velocity structure of the GRB progenitor star at the time of explosion. The host galaxy is a 
prolific star-forming galaxy with a SFR of ~40 Moyr"'. 

Key words. (Stars:) Gamma-ray bust - Techniques: spectroscopic - Stars: Wolf-Rayet - Galaxies: starburst - Cosmology: observations 



1. Introduction 



Send offprint requests to: A.J. Castro-Tirado, e-mail: ajctOiaa. es 

* Based on observations taken with the ESO's 8.2m Very Large The afterglows of long-duration Gamma-Ray Bursts (GRBs), 
Telescope in Chile. which are linked with the explosions of massive stars (see 
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Wooslev & Bloonil 2006 . for a recent review), are the most lu- 
minous optical sources in the universe for short periods of time 
( Kann et al. . 2007 : Bloom et al. . 20091) . Low-resolution optical 
spectroscopy was initially only usable to determine the redshift 
and th us place them at cosmological distances dMetzger et al.l 
[mi). Deeper insight came with the first medium-resolution 
spectrum, obtained with Keck ESI of the afterglow of GRB 
000926 (ICastro et al.L l2003h . The first true high -resolution 



echelle spectra were obtained for GRB 020813 ( Fiore et al. 



2005|), but they were of low signal-to-noise ratio. The first 



echelle spectra with good S/N were fin ally obtained for G RB 



021004, the focus of this work (see also Fiore et alil2005h . 



Such spectroscopy allows de ep ins ight into the environ- 
ments of GRBs (Prochaska et al.l l2006l) . Some highlights in- 
clude the possible detecti on of a G alactic superwind in the host 
galaxy of GRB 030329 jThone et al., 2007), and variable ab- 
sorption lines which result from direct UV pumping by the 



luminous GRB afterglow ('Pessauges-Z ayadskv et al. . 2006t 



[Vreeswijket al., 2007; D'EHaetal., 2009^). These detections 
have recently been possible for more afterglows due to the rapid 



locali zation capabilities of the Swift satellite dGehrels et al 



2004ft in combination with the rapid-response mode (RRM) 
which is now available for the Ultraviolet- Visual Echelle 
Spectrograph (U VES) at the Very La rge Telescope (VLT) (e.g., 
Vreeswiik etdl 12007: D'El ia etal.l l2009). Recently, the cov- 



ered wavelength region has been expanded all the way from the 
ultraviolet into the K band near-infrared (nI R) by the second- 



gener ation instrument X-Shootei^t VLT (iD'Odorico et al 



20061 see lde Ugarte Postigo et al.ll2009l for a first result) 



GRB 021004 was detected at 12:06: 14 universal time (UT) 
on 2002, October 4 with the gamma-ray instrument FREGATE, 
the wide-field X-ray Monitor (WXM) and the soft X-ray cam- 
era (SXC) aboard the High-Energy Transient Explorer (HETE- 
2) (Shirasakietal., 2002). GRB 021004, was a moderately 
bright, long-duration (Tgo = 100 s) event with fluences of 
6.4 X 1 0'^ erg cm'^ (7-30 k eV) and 2.3 x 10 



keV) dBarraud et aUl2002h . 



^ erg cm-2 (30-400 



The GRB was rapidly localized in flight and the position 
was reported in less than a minute, allowing rapid ground-based 
follow-up which revealed t he presence of th e fading optical 
afterglow of GRB 021004 dFox et all l2003h . This prompted 
follow-up observations at many observator ies, wh ich led to 
an exte nsive long -term coverage at X- r ay Sako & H arrison 



(l2002a|ibl). rad io dFrail & Bergeii 120021: iBergeret all 12002 : 



Poolevl l2002a"b'Pl millimeter ("de Ugarte Postigo et al."2005'). 
near-IR (de Ugarte Postigo et al., 2005; Fynbo et al., 2005) 



and optical wavelengths, both gro i ind-based dFox etaU 120031 



Bersier et al.', '2003"; 'Uemura et al.', '2003'; 'Holland et al.', '2003 
Pandeyetal., 2003; Mirabal et al., 2003; Kawabata et al.', 
2004; de Ugarte P ostigo et al.L l2005l) as well as space-based 
(iFvnbo etal.. 12005]) ^ 

The isotropic energy release during the prompt emission of 
this G RB was modest, with log Eiso - 52.65^° dKann et al. . 
boio "). and the dust extinction for this somehow reddish af- 
terglow (R-K ~ 3) was also low as is typical for many well- 



See http://www.aoc.nrao.e du/~dfrail/grb021004.dat| for the com- 
plete VLA data set. 



observed afterglows (IKann et al.L l2006l) . Despite the low en- 
ergy promptly released, this is among the most luminous af- 



terglows ever detected dKann et aL . 20061), even in compari 



son with a much larger Swift-exa. sample ( Kann et al. . 2010l) . 



The multi wavelength temporal evolution of the GRB 021004 
after glow can be explained by multiple energy injections 
(Bi ornsson et a l.. 2004; de Ugarte Postigo et al., 200 5i), never- 



theles s other scenarios can not be discarded dUazzati et al 
2002h . GRB 021004 remains one of the most well-observed 
afterglows ever 

Early-time low and medium resolution spectroscopic ob- 
servations allowed a redshift of z = 2.33 to be de- 
termined on the basis of Ly-g absor ption and emis- 



sion lines ( Chornock & Filippenkol 2002 ). Several absorp 



tion systems with o utflows velocities of few 1000 km s 
were also reported dSalamanca etal. . |2002t I Savaglio et al 



2OO2I) and studie d in detail by ^0llere 7ai.l d2002bh 



Wang et al 



(2003), 



('2003'), 'Matheson et al." (2003), Scha efer et al 
Mirabal et al.. (.2003.). .Starhng et al.. (.20051) and 



Lazzati et al.l ( 2006l) . ljakobsson et al.l ( 2005 ) also report the de 
tection of the host galaxy Lya line in narrow-band imaging. 

Fiore et al. I (12005 !) have constrained the ionization parame 



ters of the various absorption components detected above Lyor 
(at an observer restframe of 4,050 A) in GRB 021004 , inter- 
preting them as density fluctuations, Uke iLazzati et al.l (120021) 
did. Within the context of larger sa mples, the UVES spectra of 
this GRB have also been studied bv lChen et al.l(l2007h (the gen- 
eral lack of wind signatures in high -resolution spectra of GRB 
afterglows), Prochaska et al.l (l20()8 ) (N V absorption lines to- 
ward GRB afterglows). IFox et al. ( 2008h (high-ionisation lin e 
systems towar d GRB afterglows) and lTeios et al.ld2007ll2009h : 



Vergani et al.l (12009) (study of Mg II foreground absorption 
systems). 

The original expectation, as detailed high quality spec- 
troscopy of GRB afterglows became possible, was that we 
would rapidly learn much about the GRB progenitors via the 
study of the complex absorption line systems they were ex- 
pected to exhibit due to ejection events leading up to the final 
collapse. This expectation has far from proven true. 

Despite the large number of long GRB afterglow spectra 
obtained since, GRB 021004 still stands out as the one with the 
most complex set of intrinsic (0-3,000 km/s ejection velocity) 
absorption systems. At the time it also held the place as the low- 
est detected HI column density, and it still ranks between the 
lowest seven found. The complexity orig inally suggested that 



it indeed represented a display of eject a dM0ller et al.L 12002a : 



Mirabal et all 120031: iFiore et all l2()05b b ut this interpretation 
was later disputed by Chen et aLr ( 2007 ). A final interpreta- 
tion of the complex systems has not yet been agreed on, and 
this object still stands our as the best candidate of a seemingly 
rare case where the signatures of events prior to the collapse 
are displayed. The rarity alone would warrant a more detailed 
discussion of the spectral features. In addition we are adding 
UVES data below 4050 ang for a more detailed discussion 
of the absorption systems (para 3.1) and near-IR VLT/ISAAC 
spectroscopy allowing a more accurate determination of host 
properties and redshift (para 3.2). We discuss our results in the 
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Observed wavelength (A) 



Fig. 1. Overall view of GRB 021004 optical afterglow spectrum. These VLT/UVES data were obtained ~0.6 days after the GRB 
and show the Lya emission line arising from the host galaxy redshift plus some of the most prominent absorption lines systems 
at (or very close to) the host galaxy redshift: low ionization lines (like Fe II, Si II) and high ionization lines (like Si IV, C IV, N 
V), all labelled in red colour. Also shown for completeness are the absorption systems for the foreground systems at z = 1 .6020 
(labelled in green) and z - 1.3820 (labelled in blue). For clarity, the original data have been smoothed. The telluric lines in the 
A-band (7600-7630 A) are noticeable. The dotted line is the error spectrum. 



light of progenitor models for GRBs in §|4l and summarize the 
work in § |5] 

For a Hubble constant of Hq - 72 km s"' Mpc"', a matter 
density Q,„ = 0.3, and a cosmological constant Qa - 0.7, the 
luminosity distance to the host is = 18.22 Gpc, and the look- 
back time is 10.42 Gyr. All errors are given at a Icr level of 
confidence for a parameter of interest unless stated otherwise. 



2. Observations and data reduction 

2.1. Optical observations 

Observations were conducted with 8.2-m Very Large Telescope 
Units 1 (VLT/UTl, Antu) and 2 (VLT/UT2, Kueyen) 
at the European Southern Observatory (ESO) in Cerro 
Paranal (Chile). We obtained optical spectroscopy of the 
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GRB afterglow starting 0.6 days after the burst using the 
Ultraviolet- Visual Echelle Spectrograph (UVES) on UT2. 
Contemporaneous nIR observations were taken using the 
Infrared Spectrometer And Array Camera (ISAAC) on UTl 
(see Table[T]for a log of the observations). 

The optical data reduction was performed using the UVES 
context running under MIDAsH. The UVES context is struc- 
tured into reduction recipes allowing for data reduction in a 
semi-automatic manner. A master bias frame was used for the 
bias subtraction. The order tracing was performed taking ad- 
vantage of a first guess cahbration solution based on the UVES 
physical model. In all cases, except for the red spectrum cen- 
tered at 8,600 A which is affected by fringing effects, the flat 
field correction was applied after the order extraction. A Th- 
Ar lamp was used for wavelength calibration. The spectra were 
extracted using an optimal extraction method, which provides 
the object signal and the variance. After the extraction the spec- 
tra were resampled to constant wavelength bins and the orders 
merged into a single spectrum. Finally, both the standard helio- 
centric and vacuum wavelength corrections were applied. 

2.2. Near-infrared observations 

For the nIR observations we used the short wavelength mode 
in low resolution, which yields a spectral coverage of 1.42-1.83 
fim and 1.84-2.56 yum in the H and K bands, respectively. The 
pixel scale of the ESO-Hawaii detector is Of.' 146/pixel. For each 
band, 50 individual exposures of 60 seconds each were taken, 
shifting the target along the slit every 5 exposures in H (3 in K) 
in order to be able to remove the sky background, a standard 
technique in infrared astronomy. We also collected dark and 
flat field images, which were used in the reduction process per- 
formed with IRAF0. After combining the groups of exposures 
and extracting the individual spectra, we calibrated them using 
the air-glow OH emission lines present in the 2 dimensional 
images. This allows an accurate wavelength calibration with 
dispersion of 4 and 8 A/pix for H and K respectively. Finally, 
the individual spectra were median-combined into a single one. 
In the case of the H band spectrum we used a G2 spectral type 
star to remove the telluric absorption bands. The spectra were 
flux calibrated by means of the JHK-band observations that we 
obtained bracketing the spectral observations. 

3. Results 

3.1. The blueshifted absorption line systems 

The VLT/UVES data set (Fig. [T]| shows a large number of 
absorption lines, as well as Lyg in emission (see Table [3] 
for a list, and also lM0ller et al] l2002bt IStarUng et all 1200^ 



^ MIDAS (Munich Imaging Data Analysis System) is devel- 
oped and maintained by the European Southern Observatory (ESO). 
www.eso.org/projects/esomidas/midas-distrib.html 

IRAF is the Image Reduction and Analysis Facility, a general pur- 
pose software system for the reduction and analysis of astronomical 
data. IRAF is written and supported by the IRAF programming group 
at the National Optical Astronomy Observatories (NOAO) in Tucson, 
Arizona, USA. 



Table 1. Journal of the VLT GRB 021004 optical/nIR spec- 
troscopic observations. All observations were obtained on 
October 5 (UT). 



UT Time 
(start) 


Instr. 


Exp. 
time (s) 


Spectral 
Range (A) 


S/N 
ratio" 


Resolution 
A/A(A) 


03:02 


UVES 


3600 


3290-4520 


1 


~ 56000 








4620-5595 


2.2 


-51500 








5675-6645 


2.7 


-51500 


03:02 


ISAAC 


3600 


14200-18300 


5.2 


- 500 


04:10 


UVES 


3600 


3050-3870 


1 


-51300 








4780-5755 


3.7 


- 50500 








5835-6805 


5.4 


- 50500 


05:15 


UVES 


3600 


3750-4650 


2.4 


-51500 








6705-8520 


4.1 


- 52300 








8665-10400 


2.2 


- 48700 


05:27 


ISAAC 


3000 


18400-25800 


2.3 


-450 



" per resolution element. 



Table 2. Individual high-velocity metal-lines systems in GRB 
021004. 



Line ID 


C2 


CI 


D 


E2 


El 


Si 111^1,206 




y 


y 


y 


y 


HILyail,216 


y 


y 


y 


y 


y 


N Vil,238-1,242 


y 








y 


Si Hi 1,260- 1,264 




y 


? 




y 


C II /1 1,334 


y 


y 








Si IV /1 1,394- 1,403 


? 


7 


7 


? 


y 


CIV ^1,548-1,551 


y 


y 


y 


y 


y 


Al II ^1,671 




y 




y 


y 


Aim ^1,855-1,863 










y 


Fe II ^2,344-2,382 










y 


z 


2.29671 


2.29935 


2.3216 


2.3275 


2.32891 


V (km/s^ 


-3,050_^ 


-2,8 Ki;^ 


-795_^ 


-262^ 


-134^ 



" y = yes; - = no; ? = doubtful, low significance due to low S/N 
ratio in the S IV region. 

* relative to the systemic velocity of the host galaxy at z = 2.3304, 
which we set as zero. 

Embedded as CI in the wide C complex (2,730-3,250 km s"'). 

'' Two subcomponents at -2,837 and -2,803 km s"' (low ioniza- 
tion) and -2,860 and -2,810 km s"' (high ionization). 
Wide complex spanning 585-1,005 km s"'. 

^ Embedded as El in the wide E complex (40-390 km s"'). 

^ Two subcomponents at - 1 50 and - 1 1 4 km s" ' . 



Jako bsson et alj|2005l on the Lya emission line). Some of the 
absorption lines are due to two for eground systems at redshifts 



z= 1.3820 and z 
see 



1 .6020 (Mille r etalll2002hl; 



20031 see iTeios et aP 120071 120091: IVergani et al 



Mirabal et ak 
20091 for fur- 



ther discussions on the foreground absorbers), but the remain- 
ing lines are from multiple systems found in the redshift range 
2.2967 <z < 2.3291. 

In Fig.|2]we show the three wide absorption line complexes 
in the redshift range 2.2967 < z < 2.3291. The absorption pro- 
files are plotted in velocity space and m arked C, D, and E fol - 
lowing the naming convention used bv lM0ller et al.l (l2002ah . 
The higher resolution of the UVES spectrum allowed us to 
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identify four narrow line systems embedded inside the velocity 
structure of the broad systems C and E. We name those systems 
CI, C2, El and E2. In the following we discuss all of those 
systems in more detail in order of increasing velocity relative 
to the host redshift (z = 2.3304, see §1X2]). 

E complex: The Lyo- absorption width translates to veloc- 
ities in the range 40^70 km s"^ while in C IV (where we 
observe the narrowest trough), the range is 40-390 km s"'. 
The difference in the high-velocity cut-offs (390 vs. 470 km 
s~') is most easily understood as a column density effect, i.e., 
the column density of the high velocity end is too low for C 
IV to be detectable. Alternatively it could be caused by lower 
metallicity in the higher velocity gas of system E, or a different 
ionization. The simplest interpretation of the E complex is that 
it is caused by gas, spanning velocities from 40-470 km s"' 
and densities decreasing from the slowest to the fastest moving 
part. Similar profiles have been observed in QSO broad absorp- 
tion line (BAL) systems, where velocities up to 0.1c have 
been found. BAL systems are interpreted as velocity broaden- 
ing due to ejection and/or radiative acceleration. Such BAL- 
like absorption troughs are not meaningfully fitted by standard 
(Voigt) profiles that assume a Gaussian velocity field. 

Embedded in the E complex are two narrow line systems at 
134 km s"' (El) and 262 km s"' (E2). The strongest compo- 
nent. El, splits equally into two subcomponents at 1 14 km s"' 
and 150 km s"' (see Fig 3) and exhibits absorption from a large 
number of both high and low ionic species (N V, C IV, Fe II, Al 
II, Al III, Si II, Si III, Si IV). The metal Une absorption of El has 
similarities to so-called damped Lyo' (DLA) systems, which in 
a few cases have been identified with high redshift, star form- 
ing galaxies (lM0ller et al.. 2002b) but it has too low HI column 
density in this sightline to be considered a DLA. Still, El is the 
system most likely to be identified with the interstellar medium 
(ISM) in the host galaxy of the GRB, whereby the velocity off- 
set between the host and the ISM (~ 134 km s"') represents the 
velocity of the local absorbing cloud within the host. High and 
low ionic species are also seen in the second system (E2), but 
at lower significance. The two rightmost dotted vertical lines in 
Fig 2 mark the systems El and E2. 

D Complex: The D complex consists of a single, wide BAL- 
like component seen in Lyo-, C IV, and Si III. Both the C IV 
(1,548 A) and the Lyff troughs give a consistent relative veloc- 
ity range from 585 to 1,005 km s"' . We find residual flux at the 
bottom of both the Lya and C IV ( 1 ,548 A) troughs indicating 
that the gas is optically thin. We also see significant variations 
in the optical thickness as a function of relative velocity; in 
particular one will notice the two sharp "dips" marking shells 
of higher column density at the leading and the trailing edge. 
Those dips are seen in Lya, in both C IV lines and in Si III 
(marked by dotted lines in Fig 2). The Si III trough is opti- 
cally thick between those two edges, but then appears to have 
an optically thin extension towards lower ejection velocities. 
However, this trough overlaps with Lya of system C, and the 
extension is found to belong to system C (see below). . 

It is remarkable that the D and E systems are separated in 
velocity by exactly the amount required to shift the C IV 1,550 
A line complex right on top of the C IV 1 ,548 A o ne in the spec- 
trum (lM0Uer et all l2002bl: ISavaglio et all |2002|) . Such a shift 



is known as absorption-absorption "line-locking" and is not un- 
common i n QSO spectra where radiative acceleration is impor- 
tant (e.g., I Srianand et al. 2002 1). The line-locking between the 
D and E systems strongly suggests that a process similar to that 
seen in intrinsic QSO absorption systems is indeed involved 
here, i.e. radiative acceleration. 



C2 CI 



E2E1 



6 - 



5 - 



4 - 



3 - 



1 - 




-3000 -2000 -1000 

Ejection welocity Ckm/s) 

Fig. 2. In this figure we have converted the observed wave- 
lengths into velocities relative to the host redshift (z - 2.3304) 
for each ionic species significantly detected. Lya is plotted at 
the bottom. To optimize the signal-to-noise ratio two lines of 
N V (1,238 and 1,242 A) and four lines of Si II (1,260, 1,264, 
1,533 and 1,816 A) have been coadded. The rest of the lines 
are all single. The absorption components around z ~ 2.3 ob- 
served at velocities up to 3,250 km s"' indicate the presence of 
fast winds, ejected by a hot and massive stellar progenitor The 
D and E complexes both contain velocity broadened, BAL-like 
components. Within the E complex there are narrow lines. El 
and E2, presumably from the ISM of the host galaxy. There are 
also narrow shells in the high-velocity system C, CI and C2. 
The dotted lines mark the positions of El, E2, the "dips" at the 
leading and the trailing edges of D, and the CI and C2 systems. 
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C Complex: As discussed above, the wide Lya trough of 
system C is partly due to Si III absorption from system D, but 
also contains Lya absorption of a wide component spanning 
the velocity range 2,730-3,250 km s"' as we detect C II, C IV, 
Si II, Si III, Al II and N V absorption in this velocity range 
(Fig.|2]i. The finely tuned velocity offset which causes a precise 
overlap between the Lyo- (C complex) and Si III (D complex) 
lines could be due to chance coincidence (although with small 
probability, see below). Given the evidence for line-locking be- 
tween the two C IV lines between systems D and E, the overlap 
is most likely the result of radiative acceleration leading to line- 
locking of two expanding shells corresponding to the C and D 
complexes. 

Embedded in C we identify two narrow systems at 3,050 
km s"' (C2) and at ~ 2,810 km s"' (CI). C2 is a high ionisa- 
tion system displaying strong N V and C IV absorption, but 
no absorption from singly ionized ions. CI is a low ioniza- 
tion system with moderate C IV absorption and clear detec- 
tions of the singly ionized species Al II and Si II, which usu- 
ally identifies a cloud optically thick at the Lyman limit (see 
TableO. Most likely C2, although farthest away from the host 
redshift in velocity space,corresponds to the part of the wind 
that is nearest to the progenitor site in physical space. This 
ionization-velocity correlation is also observed in QSOs: the 
high ionization lines are often seen at the highest relative ve- 
lociti es for the Zabs ~ Ze m absorption systems in QSO spectra 
fe.p.. lM0ller et al.l[l994 l). As for El, CI splits into two sub- 
components, detected at 2,803 km s ' and 2,837 km s ' (low 
ionization) and at 2,810 km s ' and 2,860 km s ' (high ion- 
ization). Again, we here have the ionization-velocity correla- 
tion seen already between CI and C2 (Fig. |2]l. Similarly to 
Fiore et al.l (l2005h we do consider CI as part of the ejected sys- 
tems simply on the basis of the low probability for it to be ex- 
trinsic, although we cannot ex clude its associati on to the neigh- 
bour galaxy detected by HST dChen et al.ll2007h . In that case, it 
should be realize d that the comparat ively low H I column den- 
sity reported by Fvnbo et al.l (l2005b would be the sum of the 
H I column of the two galaxies, such that they both have even 
lower H I columns. 



3.2. The near-IR emission lines 



Assuming a redshift of the host of z = 2.3351 ("Meille r et al 
l2002 b. based on Lya emission) we searched the ISAAC spectra 
for redshifted lines of H I (Ha, and [O III] (4,959 A/5,007 
A). There are tentative detections at low significance of Ho' and 
the [O III] lines. In order to confir m the reality of the lines we 
used the method we introduced in ljensen et alj ( 200 lb . Briefly 
the method involves rebinning the spectra into redshifts for the 
rest wavelength of each expected line, in our case the lines of 
Ha and [O III] (4,959 A/5,007 A). The three "redshift spectra" 
were then coadded (with weights 2, 1 and 1 because of the 
slightly higher S/N of the tentative Ha detection). The resulting 
combined line is clearly detected, as shown in Fig. ID 

To this combined line we fit a Gaussian (overplotted as a 
blue dotted line) and we measure the systemic redshift to be 
z = 2.3304 + 0.0005, which is significantly dififerent from the 



values derived from the Lya emission line (which may be ab- 
sorbed or modified by resonant scattering) or from the absorp- 
tion systems (which are affected by the local velocity of the 
absorbing cloud in the host galaxy). The width of the Gaussian 
is found to be exactly the spectral resolution, so the lines are 
unresolved and we can place an upper limit on the intrinsic line 
width of 600 km s ' FWHM. 



J 



-1 — I — I — r- 



-I — I — r 

CI 




_1_ 



-2900 -2B00 -2700 

Ejection i^locity (km/s) 
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-300 -200 -100 
Ejection velocity (Km/s) 

Fig. 3. A zoom in of the high-resolution VLT/UVES spectrum 
at the location of the CI and El narrow-line absorption sys- 
tems. The data presented here are the co-add of the different 
ions in velocity space. Thus, CI has two subcomponents at 
2,803 km s"' and 2,837 km s"' (low ionization. Si II + Al II, 
top) and 2,810 km s"' and 2,860 km s"' (high ionization. Si III 
+ C IV, bottom ). The io nization-velocity correlat ion also seen 



in QS Os (M0ller et all 11994) and WR nebulae dSmith et al 



19841) is noticeable. El equally splits into two subsystems at 
1 14 km s"^ and 150 km s The data shown are the co-add of 
different ions (N V -i- Si II -H Al II -i- Al III) in velocity space. 
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4. Discussion 

Following the suggestion that there is often a link be- 

tween long-duration GRBs and core -collapse supernovae 




21500 22000 

Uawe length (Angstrom) 




16000 1S500 17000 

Uawe length C Angstrom) 

Fig. 4. Bottom panels: A blow-up of the GRB 021004 nIR af- 
terglow spectrum. These VLT/ISAAC data were obtained ~0.6 
days after the GRB and show the Ha and [O III] emission hnes 
arising from the host galaxy. Top panel: the co-addition of the 
three redshift spectra (as explained in the text) is shown, pro- 
viding an independent measurement of the systemic redshift z 
- 2.3304 based on the nIR line emission. In the middle and 
bottom panels we show the optimal fits (black), and the data 
itself (blue dotted). 



Using a redshift of 2.3304 and the width of the resolution 
profile, we now fitted Gaussian profiles to the data (with con- 
tinuum level and amplitude of the Gaussian as the only free pa- 
rameters) using min-square deviation . For the two [O III lines 



we further imposed a 1:3 flux ratio (lOsterbrock 



(ll989h Y 

found the following fluxes: flux = (4 + 1) x 10"^^ erg cm 



we 

-2 



S-' ([O III] 4,959 A), (12 = 
5,007 A) and (11 + 2) x 10 
is consistent with zero. 



: 3) X 

erg cm 



10"'^ erg cm-2 s"' ([O III] 



s ' (Ha). The H/? flux 



2001 iHjorfli et all l2003h as proposed in 1993 



(IWooslev-, 19931) we know that GRBs arise as the end stages of 



massive-star evolution. As we describe below, the results pre- 
sented here provide independent evidence for a massive star (> 
40 Mq) progenitor and give additional insight into the imme- 
diate surroundings and prior evolution of the star leading up to 
the catastrophic event. 

4.1. A LBV / Wolf-Rayet progenitor 

The probability that the three absorption systems CI, C2, and 
D lie in front of the GRB host galaxy as a chance coincidence 
is very small and can be computed using the number density of 
intervening C IV absorbers at 1.8 < z < 3.5 of 1.5 + 0.2 per 



unit redshift jSargent et al.Lll998h . We find that the probability 
to find three of them within the reds hift range A z - 0.0337 



Fiore et al.l (120051) used 



spanned by CI, C2, and D is ~ 10~^. 
a similar argument but added t he column densitie s to obtain an 
even smaller probablility. Like iFi ore et al.l (l2005h we conclude 
that the absorbers must be intrinsic to the host galaxy. The large 
blue shift then makes it unlikely that the physical location can 
be anywhere else than close to the GRB progenitor. 

Furthermore, the probability that the C IV 1,550 A line for 
complex E falls on top of the C IV 1,548 A line for complex 
D and the Lya line for system C falls on top of the Si III line 
for system D by chance is small, of order a few percent for 
each, although difficult to quantify due to the complex profiles. 
Therefore, it makes sense to seek a physical explanation for the 
complex structure of the C, D, and E systems. 

The multiple absorption line systems found in the spec- 
trum of GRB 021004 are naturally explained by multi- 
ple shell structures formed by the stellar winds of a mas- 



sive progenitor star (iSchaefer et all 12003). Such line-driven 



winds are expected for ver y massive stars (i Castor et al.L 11975 



Garcia Segura et al.L 119961) t hat eri d their lives as Wolf-Rayet 



stars (WR) (Ide Koter et al.L 119971) . after passin g through an 



unstable Luminous Blue Variable (LBV) phase oLanger et al 



19991) . Stellar evolution calculations for massive stars predict a 



number of different phases in which the wind properties associ- 
ated with those phases, i.e., the velocities and mass- l oss rates, 
become strongly time-varying dLanger et al. , Il999t iMaederl 



19831). The minimum initial mass for a star to become a LBV 



IS ~ 40 Mq dLanger et all 1 1994 iMaeder and Mevnetl l2000l) . 
During the course of their evolution, starting as O-class stars (O 
LBV — » WR — > SN), stars have distinct winds in which their 
velocities are proportional to the escape velocities. Thus, for 
each transition where a fast wind (i.e., when the stellar radius 
was small) follows a slow one (larger stellar radius), a shell of 
swept-up material is formed. Within this simple scheme, com- 
putations of the circumstellar medium around a 60 Mq s tar at 



the zero-age main sequence cGarcfa Segura et al.L Il996l) may 



produce the number of distinct compressed shells in the kine- 
matic range seen in Fig.|2l 



8 



A. J. Castro-Tirado et al.: GRB 021004: Tomography of a gamma-ray burst progenitor and its host galaxy 



GRB 

„ O . 
Progenitor 




Fig. 5. A sketch showing the absorption line complexes C, D and E in the line of sight to GRB 021004, not drawn to scale. 
Velocities are given in km s"' . Complex C, with the highest relative velocity, is most likely nearest to the progenitor site in space 
and probably formed at the WR phase. The high-ionization sub-component C2 is the closest one to the progenitor whereas the 
colde r sub-component C 1 is further out and we cannot exclude that it could be even related to the neighbour galaxy detected by 
HST dChen et al.1 l2007l) . Complex E has the lowest relative velocity, and the broad part of the E complex is most likely caused 
by the oldest wind formed when the progenitor evolved from being an O star into the LBV phase or in the early LBV phase. The 
narrow El component has the properties expected for the ISM of the host galaxy. We believe that the D complex is caused by 
winds formed during the unstable LBV phase. 



During the LBV phase, when hydrogen shell burning oc- 
curs, the stellar models are affected by several instabilities that 
cannot yet be modelled self-consistently and thus the unsta - 



in a number of galactic ring nebulae surrounding Wolf-Rayet 



ble LBV phase can not be fully resolved ( Langer et al. . 19941) . 



However, several LBV observations show that a number of dif- 
ferent swept-up shells can be formed. In particular, the D com- 
plex discussed above is akin to the Homunculus nebula around 
T] Car, which cont ains ~ 1 Mp and expands in the range 700 and 
1,000 km s ' ( Humphrevs and Davidson. 1994). As the forma- 
tion of each of these swep t-up shells occurs w hen the LBV 



reaches the Eddington limit (iLanger et al. [ Il999h . massive stars 



are able to form several compressed shells before ending their 
lives. 

We note that, since the LBV stellar wind is asymmetric 
and often leads to a bipolar or at least asymmetric nebula 
dOruendl et al. [ |2000h . the GRB has to be beamed in the same 
direction as the wind for us to see the shells in absorption. 
Similar, but less spectac ular, complex abso rption h as been de- 
tected in GRB 030226 dKIose et all [2 004, but see IShin et al 



2006h and GRB 050505 dBergereLal., 2006); this effect is 



not s een, however, in se veral other cases with similar quality 
data dChen et al.L l2007h . It is possible, therefore, that GRBs 



where no such complex velocity structure is seen have a dif- 
ferent orientation of the wind and GRB jet. Highly ionized gas 
(~ r = 10^ K) like CIV, N V or Si IV, similar to that found 
in the environment of GRB 021004, has already been reported 



stars (iBoronson et alilI997n . To form such gas shells, forward 
shocks with ~ 10^ kms~\ or reverse shocks with ~ 10^-10^ 
kms ' are required. The most natural explanation for the C 
complex is that it has been formed at the WR phase, with CI 
and C2 lying in an expanding free wind, before reaching the 
shocked stellar wind. In fact, the ionization- velocity correlation 
seen in the systems CI and C2 has also been det ected in WR 
nebul ae, but at lower velocities (~ 100-150 km s ' ; 
1984 . 



Smith et al 



In addition to the detailed picture of the spatial and velocity 
structure of the GRB progenitor star at the time of explosion 
(Fig. |5]l our observations have also provided insight into the 
physical mechanism producing the progenitor shell structure. 
The fact that the CI and C2 profiles are not P-Cygni, but almost 
symmetric profiles, rules out that they are very close to the burst 
(within ~ 0.2 pc) and being accelerated by the GRB ionizing 
flux to the observed velocities, as has been proposed elsewhere 



(Schaeferet al.. 2003: Mirabal et al. . 2003). The detected line 



locking also favours radiative acceleration by the progenitor 
star. 



4.2. A starburst host galaxy 

An extre mely blue hos t galaxy has been revealed in la t e time 
imaging dFvnbo et al.l 120051: Ide Ugarte Postigo et aI.L 12005 
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Jakobsson et al.Ll2005l) . In order to derive the star-forming rate 



(SFR), we have studied the emission lines seen in the combined 
optical-nIR spectrum. 

The emission line parameters are measured with a Gaussian 
fit to the emission line and a flat fit to the continuum. If all 
the He emission is attributed to star formation in the host 
galaxy we can compute the SFR as SFR (Moyr"') = 7.9 x 
10 (grg/s) (iKennicutetallil994 . In our assumed cos- 

mology, the measured Ho- intensity transforms into a Ha lumi- 
nosity Lna = (4.9 + 0.9) x lO'*^ erg s^'. This implies a SFR 
(without any corrections) of 40 + 7 Moyr"', which is much 
larger than the present- day rate in our Ga laxy. From the total 
measured Lya intensity (ICastro-Tirado et a l.. 2002), we derive 
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i42 



erg_ 



which im- 



a Lya luminosity L^ya < (9 ± 1) x 

plies a SFR < lOMpyr'^ in agreement with lFvnbo et alj (|2005[) 



and Jakobsson et al.l (l2005h . The Lya/Ha ratio of 2 implies a 
low dust content in the GRB host galaxy, as also in ferred by 
Ide Ugarte Postieo et all (l2005h and lKann et"aD (l2006 h from the 
afterglow SED; there is also no excess absorption beyond the 
(low) Galactic value dete cted in Chandra X-ray observations 
( Sako & Harrison'. '200 2a|). This is a large value for the SFR 
(see also Castro-Tirado et al.l (|2007|)), although the signal-to- 
noise ratio of the nIR spectrum implies a large uncertainty. The 
derived SFR implies that most of the ongoing star formation is 
unobscured, whic h is also in agree ment with the SCUBA re- 



sults at 850 jum ( Tanvir et all 120041) . 



The Lyes' photons produced in the ionized nebula suffer a 
continuous resonant scattering in the presence of HI. In a dust 
free nebula this mechanism effectively traps the Lyes' photons 
producing a typical P Cygni line profile in emission and ab- 
sorption, depending on the details of the density and velocity 
field. However, if even a little dust is present, most Lya photons 
will disappear by heating dust grains, because radiation absorp- 
tion by dust has a maximum cross section near the wavelength 
of Lya. Thus, the observed complex Lya profile is the result of 
both mechanisms, resonance scattering which sets the profile 
shape, and dust absorption which further decreases the inten- 
sity of the emission line. The Lya P Cygni profile can be under- 
stood in terms of these mechanisms acting during the evolution 
of the swept up HI supershell produced by a massive starburst 
in the host galaxy of GRB 021004. It co rresponds to stage 4 



in the evol utionary scheme develop ed by iTenorio-Tagle et al 



19991) and lMas-Hesse et al.l (|2003|); an stage compatible with 



GRB 021004 being the evolution endpoint of a WR star. 

5. Conclusions 

We find several absorption line groups spanning a range of 
about 3,000 km s ' in velocity relative to the redshift of the 
host galaxy. The absorption profiles are very complex with both 
velocity -broadened components extending over several 100 km 
s"' and narrow lines with velocity widths of only ~ 20 km s"' . 
By analogy with QSO absorption line studies, the relative ve- 
locities, widths, and degrees of ionization of the lines ("line- 
locking", "ionization-velocity correlation") show that the pro- 
genitor had both an extremely strong radiation field and several 
distinct mass loss phases (winds). These results are consistent 
with GRB progenitors being massive stars, such as LBVs or 



Wolf-Rayet stars, and provide further insight into the nature of 
these progenitors and their immediate environments. 

The host galaxy is a prolifically star-forming galaxy at a 
systemic reds hift z - 2.33 04, with a SFR of -40 Mpyr^' as 



also fo und bv lFvnbo et alj (l2005h and Ide Ugarte Postigo et al 



(I2OO5I), reinforcing the potential association of some GRB with 



starburst galaxies (Christensen et al. (2004): iGorosabel et al 
(l2005h and references there in). 

The ^wj/f mission with a predicted lifetime of ten years 



dGehrels et al.L 2004) will certainly bring us the opportunity 
to caiTy out high-resolution spectroscopy for dozens of future 
GRBs, and to set physical/chemical properties common to all 
GRB outflows. 
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List of Objects 



Table 3. Line Identifications in the GRB 021004 OA spectra. 



'GRB 021004' on page 1 
'GRB 021004' on pageH 



'IvacuHm (observed) (A) 


Line ID 


EW(Ai^ 


z 


4008 


HlLycr 1215.670 




2.2967 


4015 


Si III 1206.499 




2.3275 


4016 


Si III 1206.499 




2.3289 


4028 


C IV 1548.203 




1.6024 


4038 


HILya 1215.670 




2.3216 


4041 


HILya 1215.670 




2.3243 


4045 


HILya 1215.670 




2.3275 


4047 


HILya 1215.670 




2.3289 


4049* 


HILya 1215.670 




2.3304 


4084X1 


N V 1238.821 


0.46 ± 0.16 


2.2967 


4097.2 


N V 1242.804 


0.35 ± 0.15 


2.2967 


4123.9 


N V 1238.821 


0.90 ±0.14 


2.3289 


4137.2 


N V 1242.804 


0.66 ± 0.12 


2.3289 


4158.5 


Si II 1260.422 


0.25 ± 0.10 


2.2993 


4195.8 


Si II 1260.422 


0.71 ±0.19 


2.3289 


4210.2 


Si II 1264.738 


0.51 ±0.16 


2.3289 


4347.4 


Al II 1670.787 


0.49 ± 0.17 


1.6020 


4399.6 


CII 1334.530 


0.38 ±0.15 


2.2968 


4403.0 


CII 1334.530 


0.79 ±0.14 


2.2994 


4639.7 


Si IV 1393.754 


1.39 ±0.83 


2.3289 


4669.7 


Si IV 1402.770 


1.29 ±0.46 


2.3289 


5104.0 


C IV 1548.203 


1.77 ±0.12 


2.2967 


5107.8 


C IV 1548.203 


0.16 ±0.03 


2.2992 


5108.5 


C IV 1548.203 


0.84 ±0.10 


2.2998 


5112.5 


C IV 1550.777 


1.15 ±0.09 


2.2967 


5116.2 


C IV 1550.777 


0.06 ± 0.04 


2.2992 


5117.0 


C IV 1550.777 


0.67 ± 0.08 


2.2998 


5138.3 


C IV 1548.203 




2.3189 


5140.0 


C IV 1548.203 




2.3200 


5143.1 


C IV 1548.203 




2.3220 


5146.7 


C IV 1550.777 




2.3188 


5146.7 


C IV 1548.203 




2.3243 


5148.5 


C IV 1550.777 




2.3200 


5148.5 


C IV 1548.203 




2.3255 


5151.7 


C IV 1550.777 




2.3220 


5152.9 


C IV 1548.203 




2.3275+2.3289 


5155.2 


C IV 1550.777 




2.3243 


5157.6 


C IV 1550.777 




2.3255 


5161.5 


C IV 1550.777 




2.3275+2.3289 


5235.4 


Si I 1572.718 


0.80 ±0.10 


2.3289 


5559.5 


Al II 1670.787 


0.23 ± 0.07 


2.3275 


5561.9 


Al II 1670.787 


1.73 ±0.12 


2.3289 


5578.0 


Fe II 2344.214 


1.01 ±0.15 


1.3796 


5651.2 


Fell 2374.461 


0.52 ± 0.06 


1.3800 


5652.4 


Fe II 2374.461 


0.16 ±0.05 


1.3805 


5652.9 


Fell 2374.461 


0.13 ± 0.02 


1.3807 


5671.2 


Fe II 2382.765 


0.86 ± 0.07 


1.3801 


5671.9 


Fell 2382.765 


0.41 ± 0.06 


1.3804 


5672.6 


Fe II 2382.765 


0.37 ± 0.06 


1.3807 


6099.6 


Fell 2344.214 


0.59 ± 0.04 


1.6020 


6101.7 


Fe II 2344.214 


0.28 ± 0.04 


1.6029 


6154.5 


Fe II 2586.650 


0.67 ± 0.08 


1.3794 


6155.0 


Fe II 2586.650 


0.50 ± 0.09 


1.3798 




Fe II 2365.552 


0.50 ± 0.09 


1.6020 


6174.2 


Aim 1854.716 


1.15 ±0.09 


2.3289 


6178.0 


Fe II 2374.461 


0.52 ± 0.08 


1.6020 


6187.0 


Fell 2600.173 


0.67 ±0.10 


1.3794 


6187.5 


Fell 2600.173 


0.44 ± 0.08 


1.3798 


6188.0 


Fell 2599.147 


0.50 ± 0.09 


1.3808 


6201.0 


Aim 1862.790 


0.73 ±0.12 


2.3289 


6202.0 


Fe II 2382.765 


0.94 ± 0.16 


1.6028 


6654.0 


Mg II 2796.352 


1.95 ±0.11* 


1.3796 


6655.5 


Mg II 2796.352 


1.95 ±0.11* 


1.3802 


6657.0 


Mg II 2796.352 


0.71 ±0.08 


1.3808 


6671.0 


Mg II 2803.531 


1.89 ±0.41* 


1.3796 


6673.0 


Mg II 2803.531 


1.89 ±0.41» 


1.3802 


6729.0 


Fe II 2586.650 


0.25 ± 0.07 


1.6016 


6730.2 


Fe II 2586.650 


0.60 ± 0.07 


1.6020 


6732.5 


Fe II 2586.650 


0.18 ±0.06 


1.6029 


6764.0 


Fe II 2600.173 


0.56 ± 0.10 


1.6015 


6765.5 


Fell 2600.173 


0.71 ±0.09* 


1.6020 




Fell 2599.147 


0.71 ±0.09* 


1.6030 


6768.0 


Fell 2600.173 


0.46 ± 0.09 


1.6029 


6788.0 


Fe II 2607.866 


0.14 ±0.05 


1.6030 


7274.2 


Mg II 2796.352 


1.11 ±0.10 


1.6015 


7275.5 


Mg II 2796.352 


0.83 ± 0.08 


1.6020 


7277.0 


Mg II 2796.352 


0.44 ± 0.07 


1.6025 


7278.5 


Mg II 2796.352 


0.83 ± 0.08 


1.6029 


7293.0 


Mg II 2803.531 


0.89 ± 0.09 


1.6015 


7294.5 


Mg II 2803.531 


0.98 ± 0.09 


1.6020 


7296.0 


Mg II 2803.531 


0.24 ± 0.06 


1.6025 


7297.0 


Mg II 2803.531 


1.00 ±0.10 


1.6029 


7803.7 


Fe II 2344.214 


0.71 ±0.28 


2.3289 


7931.8 


Fe II 2382.765 


0.46 ±0.21 


2.3288 


16520" 


[O III] 4960.295 


-(16 ±7) 


2.3304 


16679^ 


[O III] 5008.240 


-(21 ±6) 


2.3304 


21863^ 


HI 6564.610 


-(37 + 10) 


2.3304 



" Gaussian profile fitting (for non saturated lines only). Star sym- 
bols refer to complex systems. 

* Emission line. Partially absorbed, hya flux = (1.5 ± 0.2) x 10"'^ 
erg cm"^ s"'. 

^ Emission line.[0 III] flux = (4 + 3) x 10"" erg cm"^ s"'. 
Emission line.[0 III] flux = (12 + 3) x 10"" erg cm"- s"'. 
Emission line. Ha flux = (11 + 2) x 10"" erg cm"^ s"'. 



